Partial nonlinear reciprocity breaking through ultrafast dynamics in a random 

photonic medium. 
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We demonstrate that ultrafast nonlinear dynamics gives rise to reciprocity breaking in a random 
photonic medium. Reciprocity breaking is observed via the suppression of coherent backscattering, 
a manifestation of weak localization of light. The effect is observed in a pump-probe configuration 
where the pump induces an ultrafast step-change of the refractive index during the dwell time of the 
probe light in the material. The dynamical suppression of coherent backscattering is reproduced 
well by a multiple scattering Monte Carlo simulation. Ultrafast reciprocity breaking provides a 
distinct mechanism in nonlinear optical media which opens up avenues for the active manipulation 
of mesoscopic transport, random lasers, and photon localization. 
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There is a strong analogy between the transport of clas- 
sical waves such as light and sound in scattering media, 
and the mesoscopic physics of quantum waves, such as 
electrons in a solid, or matter waves in an optical speckle 
potential Whereas in many quantum systems coher- 
ence is affected by inelastic processes, interactions are 
generally weak for light and coherent effects can be stud- 
ied in large-scale systems. The role of interactions and 
decohcrcncc in many-body quantum systems can be sim- 
ulated by introducing a nonlinear optical response 043- 
Recent pioneering experiments have explored nonlinear 
light scattering in Kerr-media and cold atomic clouds in 
the stationary regime 
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In this Letter, we demonstrate partial breaking of the 
reciprocity symmetry in a random scattering medium on 
ultrafast time scales through the dynamics induced by a 
high-intensity, femtosecond pump pulse. As a probe for 
the breaking of reciprocity we use the coherent backscat- 
tering effect, a manifestation of weak localization for op- 
tical waves. Coherent backscattering (CBS) is the con- 
structive interference of reciprocal light paths, resulting 
in a cone of enhanced intensity around the backscatter- 
ing direction [l^l . The center of the cone results from the 
constructive interference of very long light paths, which 
makes CBS a sensitive probe for multiple scattering wave 
transport. 

The reciprocity breaking effect is induced through ul- 
trafast dephasing of light paths occurring in strongly 
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scattering nanostructured semiconductors 
procity breaking requires the dynamics to be fast on the 
time scale of the photon dwell time in order to achieve 
an asymmetry between the direct and reciprocal paths. 
Such a regime of adiabatic control over light has been 
achieved only recently in high-quality photonic crystal 
nanocavities US] , and is of interest for applications in 
controlled storage and release of optical information. Op- 



tical devices exploiting concepts from mesoscopic trans- 
port are currently receiving increasing interest 16- ijj. 
A degree of active control over the flow of light in these 
devices would be of interest for applications. 

The reciprocity breaking was measured in a slab of 
porous GaP of 81 ± 2 /j,m thickness, obtained by etch- 
ing of a GaP wafer. Pump-probe nonlinear CBS cones 
were measured in a beamsplitter configuration as shown 
in Fig. [TJa). The frequency-doubled output from a re- 
generatively amplified Ti:Sapphire laser (Coherent RegA, 
200 ± 10 fs , 250 kHz) at 400 nm wavelength was used 
as a pump for inducing an ultrafast nonlinear response. 
As a probe the signal output from a parametric ampli- 
fier at 630 nm wavelength was used. Both beams were 
focused to a 25-fiin diameter spot onto the sample using 
a lens with a focal length of 15 cm to achieve a sufficient 
pump flucncc for nonlinear modulation, resulting in an 
angular resolution of 20 mrad. Backscattcrcd light was 
collected using a Si-photodiodc with a 0.5-mm aperture, 
mounted on a translation stage positioned 6.5 cm away 
from the sample. Linear polarization filters were used 
to select either the polarization conserving {R//) or the 
nonconserving {R±) scattering channels. The diffuse and 
CBS intensities R± and i?//, as well as the differential 
pump-probe signals AR± and Ai?// , i.e. the difference 
of the probe signal in presence of the pump with respect 
to R± or R//, respectively, were measured using lock- in 
detection. All signals were normalized to the diffuse re- 
flectance R± to divide out the angle-dependence of the 
setup and the (Lambertian) diffuse intensity distribution. 

For the interpretation of our experimental results, 
Monte Carlo (MC) numerical simulations of the tran- 
sient response were produced by spatiotemporal track- 
ing of photon trajectories in a finite slab. Monte Carlo 
is a technique suitable for simulating diffuse transport 
and coherent backscattering of light [20[ . In the current 
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FIG. 1: (a) Experimental setup for time-resolved pump-probe 
coherent backscattering spectroscopy, (b) Principle of reci- 
procity breaking, where the symmetry of reciprocal multiple 
scattering path (arrows) is changed by a delayed pump pulse 
(t2 > ti). 



model, illustrated in Fig. [Ub), the pump pulse induces 
a change in the complex refractive index at a time t-2,. 
In all simulations in this Letter, we used one set of op- 
timized parameters. These include the a-priory known 
parameters of sample thickness, laser pulse width and 
instrumental resolution. An internal reflection proba- 
bility of 0.6, corresponding to a linear refractive index 
no = 1.47, was includedjfollowing measurements on sim- 
ilar porous GaP slabs [21|. Parameters extracted from 
comparison between numerical and experimental data 
(shown below) include the (linear) transport mean free 
path £e = 0.45 ±0.1 /im (corresponding to k^l^ = 4.5, 
with /co = 27r/630 nm~^ the wavevector of the probe 
light in vacuum) and the (linear) transport mean free 
time te = 5 ± 1 fs, the spatiotemporal refractive index 
n(z, i) = TiQ + An(0, t), where the nonlinear index change 
An(2,t2) = (1-4 X 10-2 1.4 X lO^^z) exp(-z/Loxc) is 
contained within a photo-excited surface region of width 
iexc = 0.2 \xm just after absorption of the pump pulse 
at time t^. The imaginary part of Art describes photo- 
induced absorption of the probe light, with absorption 
time Tabs = no/(2cfcoImAn) = 0.17 ps at z = and 
time ti- Subsequent spatiotemporal dynamics involve a 
rapid, picosecond expansion of the excited volume due 
to carrier migration, with a concomitant decrease in the 
nonlinear refractive index change An, such that the total 
amount of excitation in the volume remains unchanged 
by this expansion. In addition, the excitation is assumed 
to decay with a time constant of 30 ps, which, however, 
hardly affects the results presented in this Letter. 

Figure [2ja) shows time-resolved reflectivity changes 
/S.R_\_/R^ obtained in the cross-polarized, diffuse scat- 
tering channel. The signal at positive pump-probe time 
delays is characteristic for the excited semiconductor ma- 
terial and resembles our earlier measurements for bulk 
GaP (dash-dotted line) [l3|. At negative delay times, an 
additional tail is present for the scattering layer. This 
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FIG. 2: (a) Differential diffuse reflectivity AR±/R± as a func- 
tion of the time delay between pump and probe pulse. Lines 
indicate model calculations oc t~^^^ (dash), same but includ- 
ing an exponential cutoff with time-constant 10.2 ± 0.7 ps 
(dotted, blue), and Monte Carlo model (thick line). Dash- 
dotted line: signal from a homogeneous (i.e. non-disordered) 
GaP substrate, (b) Angle- and time-resolved differential re- 
flectivity AR^^/R± in the polarization conserving channel, 
for negative (red) and positive (blue) time delays ti —t2. (c) 
Same as (b) for several delay times, including results from 
MC model (red lines), (d) CBS intensity traces without pump 
('linear') and with pump reconstructed using the differential 
response curves from (c). 



ing the sample before the pump and retained inside the 
sample for thousands of scattering events before under- 
going photo-induced absorption initiated by the pump 
pulse. 

The negative-time tail of Fig.[2Ja) thus acts as a probe 
of the distribution of photon dwell times in the scattering 
medium, complementary to time-resolved interferometry 
experiments j22], [2^ . Generally, the pump-probe signal 
can be expressed as a convolution of the diffuse path- 
length distribution ad{t) and a photo-induced absorption 
term F{t), yielding 



AR^{t) 
R± 



1 

Ttot 



dhad{t2-t)F{t2), (1) 



tail corresponds to the contribution from photons enter- 



where t = ti — 12 is the pump-probe delay time. The total 
albedo aj^°* results from integrating ad{t) over time. In 
the diffusive limit, the path length distribution is given by 
ad{t) ex for times larger than the extinction mean 

free time l2J]. For the simple case of diffusion and a 
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FIG. 3: DifTerential reflectivity for parallel (open circles, red) 
and cross-polarized (closed triangles, black) channels normal- 
ized to their values eX Q = —270 mrad. Closed diamonds, 
green: Linear coherent backscattering cones. Lines: results 
from Monte Carlo model for linear and differential reflectiv- 
ity cones. For negative delay times, the differential reflectivity 
in the parallel channel is larger than 2, which indicates the 
occurence of partial reciprocity breaking. 
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FIG. 4: (a) Full width at half maximum values of the CBS 
cone width, with calculation using the CBS albedo Eq. ([2]) 
(dashed line, red) and Monte Carlo model (thick line), (b) 
Experimental pump-probe enhancement factors (circles) in 
comparison with the Monte Carlo model for the full, complex 
nonlinear refractive index as described in text (line) and using 
only the imaginary part of An (dashed line, red). 



step response F(t) — r]oQ{t), Eq. ([T]) can be solved ana- 
lytically, yielding ARj_{t)/R± = ryolV^el"^/^ for t < -U- 
This relation is shown by the dashed line in Fig. ^a). 
Clearly, this curve overestimates the contribution of long 
paths in the tail for t < —2 ps. This discrepancy can 
be corrected by including an exponential cutoff in the 
distribution ad{t) with time constant of 10.2 ± 0.7 ps 
corresponding to 2 x 10^ scattering events [dotted line 
(blue) in Fig. [UJa)]. The loss of long light paths in the 
pump-probe signal is related to the finite system size as 
explained below. Our measurements show that pump- 
probe transient absorption forms a sensitive probe of very 
long light paths in the photonic medium. 

The above form of F{t) is valid under the assumption 
that only a thin surface layer is affected by the pump 
pulse, which is a good approximation given the bulk ab- 
sorption length of GaP of around 0.2 /im at 400 nm wave- 
length |25|. Therefore, the relative amount of absorption 
is approximately constant {r]o ~ 10~^) for all probe pho- 
tons which undergo a scattering sequence whose length 
exceeds the pump-probe delay time t. This picture is con- 
firmed by the MC model (thick solid line). It was found 
from the MC model that the slope of the negative-time 
tail is determined entirely by the 25-/xm transverse spot 
size of the pump beam. Long light paths diffuse out of 
the excited region, resulting in the observed reduction of 
the long-time tail. The spatiotemporal expansion mainly 
affects the decay at positive times in Fig. ^a) . The real 
part of An does not affect the diffuse tail; however it is 
mainly responsible for the reciprocity breaking effect as 
shown below. 

We proceed to investigate the effect of ultrafast excita- 
tion on the CBS intensity in the polarization conserving 
channel. Figure [DJb) shows the experimental differential 
reflectivity AR///R±. We emphasize that these nonlin- 
ear CBS traces are obtained as the difference between 



two CBS intensity cones respectively in presence and ab- 
sence of the pump pulse, as is typical for lock-in detec- 
tion. Therefore, the differential reflectivity corresponds 
to the intensity which is removed from the linear CBS 
cone as a consequence of photo- induced effects. This can 
be observed in Fig. [2Ic,d) where we have used the differ- 
ential signals to reconstruct the modified intensity CBS 
cones at various delay times. The cone rounds off and ef- 
fectively broadens as long light paths arc removed. Both 
the coherent backscattering and the diffuse background 
show a maximum decrease of around 12% around the 
pump-probe delay of ps. 

The contribution of reciprocity breaking to the overall 
pump-probe signal can be analyzed by comparing the ra- 
tio of of the intensities removed from the CBS with those 
removed from the diffuse albedo. In particular, the more 
the total CBS cone decreases due to reciprocity break- 
ing, the more the differential CBS cone increases, which 
may thus assume a relative height of larger than 2. Fig- 
ure [3] shows the three differential CBS cones i, ii, iii of 
Fig.[2)Jc) (dots, red), where the change in the diffuse back- 
ground at -270 mrad is normalized to 1. For comparison, 
also the diffuse pump-probe signals are shown (triangles, 
black). In comparison to the linear CBS intensity cone 
(diamonds, green), the differential pump-probe cones are 
narrower and have a higher amplitude. The values for 
the angular width and differential cone height are shown 
in Fig. |DJa,b) against pump-probe delay time. 

To understand the trends in Fig.[3]and|4l we start from 
the characteristic time- and angle-dependent albedo of 
the CBS cone given by [23| 



aciO,t) ~ ad{t)e' 



(2) 



with 6 the backscattering angle. At exact backscattering, 
Eq. ^ reduces to the diffuse albedo, ac(0,i) = C(d{t). 
The angular range 9 over which paths of a characteristic 
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time duration tc contribute to the CBS-cone is given 
by d{Tc) = {Me/kltlTcf/'^. This relation shows that 
the angular width is reduced proportional to the square 
root of the path time duration. At long negative delay 
times t = ti — only long light paths with dwell times 
exceeding \t\ contribute to the nonlinear signal, and hence 
Tc = \t\. This relationship - which clearly explains the 
fact that the nonlinear CBS traces are narrower than the 
linear ones, as already observed in Fig. [3] - is indicated 
by the dashed line (red) in Fig. IDJa). The calculated 
curve includes the instrumental resolution of 20 mrad, 
and was fitted to the experimental data at negative delay 
times in order to determine tg- Around i = 0, the cone 
width exhibits a maximum, which is reproduced by the 
MC simulation (thick solid line) and can be attributed 
to the spatio-temporal dynamics of the excited surface 
region after excitation. In the MC simulation, the latter 
is modelled as an expansion of the excited surface from 
-^exc = 0.2 /im to 2.25 ^m within a time of 2 ps. Due to 
this expansion, the size of the excited region encountered 
by the probe photons increases for positive delay times, 
leading to stronger absorption of long scattering paths, 
and consequently, to a narrower differential cone. 

The maximum height of the Ai? // / Ai?j^ enhancement 
factor is shown against delay time in Fig. |4|b). If this 
factor is larger than 2, this indicates that relatively more 
intensity is removed from the CBS- than from the diffuse 
scattering-albedo. Figure|3^b) shows that this additional 
contribution is only present at negative delay times, i.e. 
for probe photons arriving before the pump pulse. This 
behavior is reproduced by our MC simulations, and can 
be attributed to the asymmetry between direct and re- 
ciprocal light paths induced by the nonlinear refractive 
index change [26|. Only changes fast with respect to the 
photon dwell time are able to break reciprocity. Such 
fast changes are mainly obtained at the leading edge of 
the ultrafast nonlinear response of GaP. As long as the 
probe photons arrive before this nonlinear edge, a suffi- 
cient reciprocity breaking is obtained to produce an effect 
in the nonlinear CBS-enhancement. Given the maximum 
value of 2.5 for the Ai?///Ai?^ ratio and our observa- 
tion that the diffuse scattering albedo of long light paths 
is suppressed by 770 — 10%, we conclude that the CBS 
albedo of these light paths is reduced, on average, by ap- 
proximately 15%, corresponding to a partial reciprocity 
breaking of about 5%. 

Dephasing and absorption both result in an asymmet- 
ric reduction of the amplitudes, and thus a reciprocity 
breaking effect. In the MC simulation (line) in Fig.|4l[b), 
good agreement with experiments is found when assum- 
ing for the real part to be ten times larger than the imag- 
inary part of An. This implies that the effect of dephas- 
ing is much stronger, consistent with our earlier results 
[13| . A purely imaginary An (dashed line) provides only 
a small reciprocity breaking effect which cannot explain 
the experimental enhancement factors. 



In conclusion, we have demonstrated partial suppres- 
sion of reciprocity in a random medium through the co- 
herent backscattering effect. Reciprocity is a fundamen- 
tal property of photonic media, and its ultrafast control 
opens up a new approach to nonlinear manipulation of 
photonic eigenstates. Our work is the first step toward 
the control of mesoscopic interference phenomena in ran- 
dom media. We achieve 5% reciprocity breaking using 
only a 200-nm thin excitation region. Typically, localiza- 
tion lengths in random media and gain lengths in random 
lasers amount to at least several micrometers [l|, [l^ • Ul- 
trafast reciprocity breaking of closed loops and random 
laser modes is therefore feasible but requires further opti- 
mization of the pumping conditions to achieve dephasing 
in a large excitation volume. This may be achieved using 
a longer pump absorption length, two-photon absorption, 
or use of two-dimensional waveguides where pump light 
can be coupled in from the top of the structure. 
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